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We calculate the pyrocoefficient, static dielectric susceptibility profiles and its thickness depen- 
dence of ferroelectric thin films. Also, the temperature dependences of above quantities have been 
calculated. For the calculations we use Landau phenomenological theory, leading to Lame equations. 
These equations subject to boundary conditions with different extrapolation length on the surfaces 
have been solved numerically. The divergency of pyroelectric coefficient and static dielectric suscep- 
tibility in the vicinity of thickness induced ferroelectric phase transition (i.e. at Z « or T ~ Td) 
has been shown to be the most prominent size effect in ferroelectric thin films. 

We also calculate the dependence of critical thickness Ic and critical temperature Td on extrap- 
olation length and film thickness. The latter defines temperature - thickness phase diagram of the 
film with paraelectric phase at / < Zc, T > or T > Td and ferroelectric phase sd, I > Ic and T < Td- 
We compare our theory with available experimental data. 
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I. INTRODUCTION 

Ferroelectric thin films attract much attention of the 
scientists due to important fundamental problems related 
to size effects resulting in drastic difference between phys- 
ical properties of the films and of the bulk ferroelectrics 
(see e.g. |^). On the other hand in the recent years a lot 
of efforts have been spared to the application of ferroelec- 
tric thin films in pyroelectric and piezoelectric sensors, 
in some types of random access memory devices, optical 
memories, microswitches etc. 1^,^. The prospect of these 
and other applications depends strongly on the solution 
of fundamental scientific problems of thin films anoma- 
lous properties and their dependence on the film thick- 
ness. The latter seems to be especially important since 
the film thickness can be easily adjusted, this adjustment 
can lead to essential change of the physical properties. In 
particular, to obtain thin film pyroelectric sensors or ca- 
pacitors with high performance one has to know the film 
thickness that gives maximal value of pyrocoefficient or 
dielectric permittivity respectively. Keeping in mind that 
recently giant dielectric response (several hundred thou- 
sands) was observed in ferroelectric thin film multilayers 
[Ql and its origin was shown to be thickness induced fer- 
roelectric phase transition one can suppose the exis- 
tence of anomalous dielectric response and other proper- 
ties peculiarities in the single film with some specially ad- 
justed thickness. The critical thickness corresponding to 
disappearance of spontaneous polarization (thickness in- 
duced ferroelectric phase transition) had been calculated 
recently in several papers (see e.g. 1^,01) in the framework 
of phenomenological Landau theory. There are only few 



works devoted to dielectric permittivity calculations (see 
e.g. ||]). To the best of our knowledge there are no calcu- 
lations of pyroelectric coefficient profile and its thickness 
dependence in ferroelectric thin films. Besides that, all 
the calculations were performed in supposition that ex- 
trapolation lengths (i.e. boundary conditions) on both 
film surfaces are the same. This seems to be rough ap- 
proximation due to different conditions on the film sur- 
faces (e.g. substrate and air). Moreover, the greatest 
part of calculations of ferroelectric thin films properties 
were performed numerically mainly for the parameters 
of BaTiOs or PbTiOs films (see e.g. |,0), which makes 
it difficult to apply the results for arbitrary ferroelectric 
film and to extract general features of size effects in the 
films. 

In the present paper we perform the analytical calcula- 
tions of pyroelectric coefficient and dielectric permittivity 
in ferroelectric thin film in the framework of the phe- 
nomenological theory. The profiles, thickness and tem- 
perature dependences of these quantities were obtained 
by analytical solution of Lame-type equations with dif- 
ferent extrapolation lengths (i.e. different boundary con- 
ditions) on the surfaces. The theory predicts the diver- 
gency of dielectric permittivity and pyrocoefficient in the 
film with critical thickness. We compare our theory with 
available experimental data. 



II. THEORY. GENERAL EQUATIONS 

In what follows we shall consider a film polarized along 
z direction normal to the surfaces of the film, i.e. polar- 
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ization P = Pz, Px = Py = 0. This type of polarization tion of an external electric field |ri| . 

can appear as a result of self-polarization of a film grown In the phenomenological theory Pz{z) = P can be ob- 

under special technological conditions without applica- tained by the minimization of free energy functional 
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where E = is an external electric field, I is the film 
thickness, Si and S2 are extrapolation lengths on the film 
surfaces. The coefficients a, (3, 7 are those of bulk ma- 
terial renormalized by internal mechanical stresses orig- 
inated from mismatch of substrate and film lattice con- 
stants, expansion coefficients, growth imperfections (see 
|p^,p|) and depolarizing field. Since this renormalization 
can result into transformation of first order phase transi- 
tion into that of a second order Eq.(|^) can be valid 
for the materials both with first and second order bulk 
phase transition. 
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FIG. 1%' The thickness dependence of Pmax and m for sev- 
eral values of extrapolation length. Curves 1, 2, 3, 4 and 
5 correspond to extrapolation lenghts (5oi, ^02) = (0.1, 0.5), 
(0.1, 2); (0.1, 10), (0.5, 10) and (2, 10) respectively. Arrows 
correspond to critical thickness. 

The minimization of the functional (|l|) implies follow- 
ing Euler-Lagrange equation 
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and boundary conditions 
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It is seen from Eqs.(^) that extrapolation lengths deter- 
mine the crossing points zi and Z2 between z axis and 
tangent lines to P{z) curve at the points z — and z = I 



respectively, i.e. 61 — —zi and 62 — Z2 — I ■ From this 
geometric definition of di and S2 one can see that (5i > 0, 
62 > when polarization on the surfaces is smaller than 
that in the film and both Si and S2 are negative in op- 
posite case. Since the conditions on two surfaces can be 
strongly different (e.g. substrates, electrodes etc.) it may 
be possible to have different signs of i5i and 62, e.g. Si > 
if P{z = 0) < P{z) and (^2 < if P{z = I) > P(z). In 
what follows we shall consider the case (5i > 0, (52 > 0. 
Eqs. (|l|) and (|2^) lead to the following equation for equi- 
librium inhomogeneous polarization in the film: 
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which has to be solved subject to the boundary condi- 
tions (2b). 

Since a = ao{T — Tc), where Tc is the temperature of 
ferroelectric phase transition of the thick film (see e.g. 
Q), inhomogeneous polarization of the film P{z) has to 
depend on T and E. This lead to inhomogeneous py- 
rocoefficient n(z) — {dP{z)/dT)p=p^ and linear dielec- 
tric susceptibility x(z) — {dP{z) / dE)p=pg where Pg is 
spontaneous polarization defined by Eq. (^ at = and 
boundary conditions (pb|) . The differentiation of Eqs. (||) 
and ( p^ ) gives following differential equations for pyro- 
coefficient and dielectric susceptibility 
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The Eqs. ( [4a| ) and (5b) define, respectively, pyroelectric 
coefficient and dielectric susceptibility profiles, i.e. 11(2:) 
and x(^)j Eis well as their dependences on temperature 
and thickness of the film. Mean values of pyroelectric 
coefficient 11 and dielectric susceptibility x can be calcu- 
lated as follows 
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n = y / Ii{z)dz, X = J I x{z)dz. 
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FIQ2 2. The dependence of the critical thickness on the 
extrapolation length. Curves 1, 2, 3 and 4 correspond to 
S02 = Soi, S02 = 0.2 (5oi, S02 = 0.1 (5oi and S02 = 0.05 (5oi 
respectively. 



III. INHOMOGENEOUS SPONTANEOUS 
POLARIZATION 

Spontaneous polarization can be found by solution of 
Eq.(0) at E = with boundary conditions (Ebh. In di- 



mensionless variables 



Pi = ^> ^1 = F' '1 = r ' '^oi'2 = (7a) 

-TO to '0 '0 



where Iq = Pq = \/-a/f3 {a = ao(T - TJ, 

T < Tc) is the homogeneous polarization of a thick film 
(where the contribution of gradient can be neglected), the 
equation for spontaneous polarization of thin film has the 
form: 
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The solution of this nonlinear equation subject to above 
boundary conditions can be written in terms of elliptic 
sine (sn(...)) (see e.g. @0) 
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where m and zo are the parameters defined by the bound- 
ary conditions and film thickness. 

In particular, value of zq can be obtained from Eqs. 
(H) and (2b) at zi — 0. It can be written as 



Zq — y/T+m F(arcsin(/i), m) (9) 
Values of parameter m can be found from the equation 



I {2K{m) — F(arcsin(/i), m) — F(arcsin(/2), m)) 
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Here 
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Note, that to obtain Eq.(g^ we used Eqs. (g) and (g^) 
ai z = l. In Eqs.(P,p^) K{m) and F{ip,m) are complete 
and inconiplete elliptic integrals of the first kind respec- 
tively ||l|,|l5|. The Eqs. (|l)-(0) are general formulas for 
calculation of inhomogeneous polarization of ferroelectric 
thin films. 

It should be noted, that nontrivial solution of equa- 
tion (^ in the temperature region T > Tc and zero ex- 
ternal field E = can not satisfy boundary conditions 
( pa| ) with positive extrapolation lengths, so only trivial 
solution Ps = is possible in this case. 



IV. THICKNESS INDUCED FERROELECTRIC 
PHASE TRANSITIONS 

The analysis of Eq. (|l^) has shown that values of m lie 
in the interval [0,1], m = 1 being the value for thick 
film. We report maximal polarization and m values 
(respectively on left and right side scales) on Fig.l as 
the functions of a film thickness. One can see that at 
I/Iq ~ SPm ~ ^Oi i-fi- maximal polarization in the film is 
close to that of a thick film. It follows from Eq. (||) , that at 
TO = polarization equals zero. This situation appears 
at some critical thickness Ic, that was found from Eqs. 
(p^), (^ at TO = 0. For the considered case Sqi > 0, 
S02 > 0. It can be written in the form 
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Eq.(y_2|) represents the dimensionless critical thickness 
at which thickness induced ferroelectric phase transition 
occurs, i.e. ferroelectric phase disappears in the films 
thinner than critical thickness Ic- 

One can see from Eq.([l^) that maximal value of he 
{he — Ti") corresponds to Sqi — S02 = 0, i.e. for 
Piz)\z=o = Pi^)\z=i = (see Eq. Increase of (5oi 

value leads to decrease he so that at 601,2 — > oohc 0. 
In the limit Sqi —f 00, S02 — > (polarization on one of 
the surfaces equals zero and acquires some maximal value 
on the other surface) he 7r/2 (see Eq.(p^). In Fig. 2 
the values of he for intermediate Sqi 7^ S02 are reported. 
In particular it follows from this Figure that Sqi ~ 40, 



$02 = 2 is far enough from the limit he ^ '"'/S, which 
can be achieved at smaller 602 values. 

The thickness induced ferroelectric phase transition 
can occur also in the film with some arbitrary thickness 
I ^ Ic &t critical temperature T = Td- Qualitatively this 
follows from the fact that m depends both on I and T 
so that m = can be obtained at some special (criti- 
cal) values of film thickness or temperature. Quantita- 
tively Tci as a function of film thickness and extrapola- 
tion length can be obtained from Eq. ( [l^ ) by substituting 

^0 = \/^/ (aoTed '- Tci)), Tci = Tci/Tc and h = I- This 
transforms Eq.(12) into that for Td, namely 



A\/l - Tci = TT - 
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are renormalized dimensionless thickness and extrap- 
olation lengths (compare ( [Za] ) and (13b)). Expression 
(13a) determines thickness dependence of dimensionless 
temperature Td of thickness induced ferroelectric phase 
transition. We show it in Fig. 3 for several values of c?i and 
d2. It is seen that the limit of thick film {Td — Tc) for all 
the curves corresponds to / > 15/0(7" = 0). However this 
limit can be achieved for thinner film { I > 2lo{T — 0)) 
in the case of large enough di and d2 (see curve 1). The 
large di,2 can be related to large 5i,2 and/or values 
and/or small 7 value (small contribution of polarization 
gradient), which can be the consequence of strong enough 
ferroelectric phase transition in a thick film. Note that 
values of critical thickness shown in Fig.l, 2 are given in 
the units Iq = y^j/ (q!o(7c — T)), while in the Fig. 3 they 
are given in the units lo{T — 0). 

Therefore ferroelectric phase transition in a film can 
be reached by variation of the film thickness at some 
fixed temperature or by changing the temperature of the 
film with fixed thickness. As a matter of fact the curves 
in Fig. 3 determine phase boundaries between paraelec- 



tric and ferroelectric phases in coordinates temperature - 
thickness of the film. Nanely, a,t I < h paraelectric phase 
(polarization P — 0) exists for all temperatures, while 
a,t I > Ic it can be both ferroelectric phase {P ^ 0) at 
T < Td and paraelectric phase at T > Td- 



V. PYROELECTRIC COEFFICIENT AND 
DIELECTRIC SUSCEPTIBILITY PROFILES AND 
THICKNESS DEPENDENCE 

A. Pyroelectric coefficient 

The solution of Eq. ( ^ ) with the boundary conditions 
( ^ ) defines the distribution of pyroelectric coefficient 11 
and has the following form (see Appendix 1 for details): 



n(.x) = Cryi{x) + CPyy2{x) + yf (x) 
where function yi{x) have the form: 
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Here IIq is a thick film pyroelectric coefficient and the 
folowing notations are introduced: 



zi + Zq 



,yo{x) = sn(a;,TO). 



In Eqs.( 14b| , 14c) E{Lp,m) and am(a;) are incomplete el- 
liptic integral of the second kind and elliptic amplitude 
function respectively ||l4,n3|. We give the explicit ex- 
pressions for Cf and Cy^n Appendix 1 due to to their 
cumberso me form. These expressions along with Eqs. 
( |l3cj - [l4a|) define the profile and temperature dependence 
of pyroelectric coefficient of a thin film. 



. 8 



. 6 



Tel 
Te 



. 4 




0.2 



10 
///o(T=0) 

FIOs'S. The thickness dependence of the critical tem- 
perature of thickness induced ferroelectric phase tran- 
sition for different dimentionless extrapolation lengths: 
(di, d2) = (10, 15), (2, 0.5), (0.2, 0,1) respectively for solid 
hues 1, 2 and 3. The dashed and dotted lines illustrate the 



approximate dependencies given by Eqs. (28b) and (28a) re- 
spectively with the same parameters as for the curve 2. 



This coefficient profiles are reported in Fig. 4 for sev- 
eral TO values. One can see from Fig.l that m can be 
considered as a characteristic of a film thickness: to w 1 
corresponds to the thickest film, while m —^ corre- 
sponds to the film with critical thickness, i.e. to the 
thinnest possible film. One can see from Fig. 4 that for 
thick enough films pyroelectric coefficient profile has two 
maxima in the vicinity of a film surfaces and a flat min- 
imum in the middle part of the film. The value of n(z) 
in this minimum depends on a film thickness: it is close 
to pyroelectric coefficient of a thick film for the thick- 
est film (see curves 6 and 7) and it increases with a film 
thickness decrease (see curve 5) up to complete disap- 
pearance of the minimum and its transformation into a 
flat maximum for some intermediate film thickness (see 
curves 3, 4). For thinner films pyroelectric coefficient 
profile has one maximum. Its height and sharpness in- 
creases with a film thickness decrease (see curves 1, 2). 
The strong increase of maximal 11(2) value for small m 
speaks in favour of possible anomaly of pyroelectric co- 
efficient in the film with critical thickness. To check this 
supposition we performed the calculations of the mean 
pyroelectric c oeffic ient 11 on the base of Eq.(^). With 
respect to Eq.(13c) one can write 11 in the following form 
(see Appendix 2 for details): 



(14d) 



where mean values of functions yi{x) have the form: 
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{x2 - xi) {1 - m) 
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(15a) 
(15b) 
(15c) 
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Here the following notations are introduced: 



trie coefficient for a film with critical thickness. 



xi = i^(arcsin(/i), m), X2 — 2K{ra) — F(arcsin(/2), m) 
ei — i?(arcsin(/i), to), 62 = 2E{m) — £'(arcsin(/2), to) 

and E{lp, to) is the complete elliptic integral of the sec- 
ond kind |lj,|l^ . For the sake of illustration we depicted 
thickness dependence of inverse mean pyrocoefficient 11 
in Fig. 5 for several (5oi, i5o2 values. One can see that at 
? = IIo/II ^ 0, i.e. there is the divergency of pyroelec- 



Temperature dependence of pyroelectric coefficient in 
the film with arbitrary thickness can be derived on the 
base of obtained expressions taking into account that 
the relation between the film thickness and parameter 
TO (see Eq.(lO)) depends on T via parameter Iq{T) — 

{ao{Tc ~ T)) = /o(t = 0)/\/l - T, T = TlTc in vari- 
ables given by ( |l3b| ). In these variables Eq.(pX|) can be 
rewritten as 



— T = \/l — m \2K{ni) — F (arcsin (/ (to, diVl — t)) , to) — 
~F (arcsin (/ (to, d2y/l — t)) , to)] 



(16) 



This expression defines the temperature dependence 
of parameter to. This makes possible to calculate the 
temperature dep endence of pyroelectric coefficient on the 
base of Eq.(14d). The results are reported in Fig. 6 by 
solid lines for several film thickness. The divergency of 
n(T) at T = Tci can be seen, the Td values for the con- 
sidered film thickness are marked by open circles in Fig. 3. 



B. Dielectric susceptibility 

Contrary to pyroelectric coefficient that exist only in 
ferroelectric phase, dielectric response is the property of 
both ferroelectric and paraelectric phase. 

Let us begin with its calculation in ferroelectric phase 
that exist at I > Ic and T < Td- In such a case it can 
be calculated similarly to pyroelectric coefficient on the 
base of Eq.(|5a,5b) (see Appendix 1 for details): 



X{x) = ClSii^) -f Ctv2{x) -f vtix) 



where 



yt{x)^-2xo{l + m) 



1 + TO — 2to yo{x) 



(17) 



(18) 



-l/2a is thick fflm susceptibility. Coefficients 



ch 



and C|'* differ from Cf^ and Cf ^ in Eq.(13c) because of 



the difference between y^^{x) and y%^{x) (see Appendix 
1). Therefore profiles of dielectric suscept ibilit y can be 
described by Eq. ([l7|) with respect to Eqs.(14a, |l|) and 
the values of Cf^ given in Appendix 1. 

The dielectric susceptibility profiles calculated on the 
base of the equation ( [l7| ) are reported in Fig. 7 for several 
values of parameter to. Comparison of Fig. 7 with Fig. 4 
shows that the profiles of pyrocoefficient and dielectric 
susceptibility looks qualitatively similar with the same 
peculiarities of x{^) for thicker and thinner films as it 
was described in previous subsection. However quantita- 
tively the behaviour of x{^) for the thinner films differs 



from that for H(z): the rate of x{^) increase with ffim 
thickness decrease is several times larger then the rate of 
H(z) increase (compare the Figs. 7 and 4). This speaks in 
favour of statement that dielectric susceptibility is more 
sensitive to size effects than pyroelectric coefficient. To 
show that the increase of x in the thinner film is re- 
lated to approaching to thickness induced ferroelectric 
phase transition we performed the calculation of mean 
dielectric susceptibility x similarly to fl calculations (see 
Eq.(ra)) 



X = t-l yi + L-2 2/2 



•ch- 



(19a) 



Integration of the third term in Eq.(p^ yields following 
expression (see Appendix 2): 



Xo 



X2 — xi \ — m 



X2-X,-^-^^^^]. (19b) 
1 — TO 



This expression along with Eqs.(15a, 15b) defines the 
mean value of dielectric susceptibility x{^)- It is shown 
on Fig. 8 by right branches (/ > Ic) of solid lines for several 
values of Jqi and 6q2 ■ It is seen that at critical thickness 
marked by arrows Xo/x = i.e. x ^ oo at Z ^ Zc so that 
we have dielectric susceptibility divergency in the film 
with critical thickness. Temperature dependence of sus- 
ceptibility can be calculated with the help of the same 
formulas (Eq.(p^) keeping in mind that parameter m 
depends also on temperature. The results of calculations 
are shown as the left branches (T < Td) of solid lines 
in Fig. 9. The divergency of x(t) at r = Td = Td/Tc 
is seen clearly, the values of Td in the considered films 
being marked by open circles in Fig. 3. Comparison of 
pyroelectric coefficient and dielectric susceptibility tem- 
perature dependencies (compare Fig. 6 and Fig. 9) shows 
that the rate of approaching to infinity as T ^ Td is 
smaller for pyroelectric coefficient than for susceptibility 
and Xo/x is linear function of T/T^ contrary to Hq/H. 

Let us proceed now to calculation of susceptibility in 
paraelectric phase, i.e. ail <lc and T — const or T > Td 
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{Tel < Tc see Fig. 3) and I = const. In such a case Ps = 
so that Eq. (|5^) gives foUowing equation for susceptibiUty 



ax ~ 7 



(20) 



which has to be solved subject to boundary conditions 

(!)• 

The equation (pQ) is ordinary differential equation of 
the second order with constant coefficients and its solu- 
tion can be expressed via trigonometric (when 7/a < 0) 
or hyperbolic functions (when j/a > 0), while in fer- 
roelectric phase the coefficient of Eq.(|5a|), including Pj, 
dependends on z and as a result its solution is expressed 
via elliptic functions rather than simple trigonometric or 
hyperbolic. 

If temperature is less then the parameter a = 



ao{T - Tc) < so that 7/a < and solution of Eq.(|2^) 
in dimensionless variables (7a) subject to boundary con- 
ditions (|^ has the form 

X{z) = - 1 1 + -r- [(cos (h) - (5o2 sin {h) -~ 1) sin (zi) - 
a { Ai 

((5o2 COS (^i) + sin (^i) -I- Sqi) cos (zi)]} . (21a) 

Here 

Ai = (<5oi + S02) cos (^i) -f (1 - S01S02) sin {h) (21b) 



The expressions (21a, pig ) define the dielectric suscep- 
tibility profile in paraelectric phase. Mean value of the 
susceptibility is found on the base of Eq. (^ and the in- 
tegration yields: 



2 (1 - cos ih)) + (Soi + 602) sin (^i) 



h {{601 + S02) cos (^i) + (1 - 601802) sin (^i)) 



- 1 



(22) 



Thickness dependence described by this expression is 
reported in Fig. 8 by left branches (/ < Ic) of solid lines 
for several values of dimensionless extrapolation lengths 
Sqi and 602- It is seen that — > Z^) — s- 00 while 
approaching 1^ both from paraelectric and ferroelectric 
phase sides, however the rate of this approaching is differ- 
ent and the value of susceptibility in paraelectric phase is 
larger than that in ferroelectric phase. Temperature de- 
pendence of mean dielectric susceptibility in paraelectric 
phase can be calculated on the base of Eq. ( p^ ) allowing 
for li — X\/l — T. The results of calculations are reported 
in Fig.9 by right branches (Tc > T > Td) of sohd fines, 



the divergency at T = Td as well as linear dependence of 
Xo/x on T /Tc being the characteristic features. 

When temperature is larger than Tc (so that a > 0, 
7/a > 0)), Zo becomes imaginary so that zi and 5oi are 
also imaginary. In this case the trigonomet ric f unctions 
(sin(. . .) and cos(. in the solution (Eq.( ^la| )) of dif- 
ferential equation (EOl) convert into hyperbolic functions 
isinh(. . .) and coshJT"!.). In dimen sionl ess variables ( |7a| ) 
the substitution — a — > a in Eqs.(pTb|,p^) gives the ex- 
pressions for distribution of susceptibility x in paraelectic 
phase T>Tc: 



J 



x{z) = — ( 1 + -r- {[cosh (/i) -I- (S02 sinh (^i) - 1] sinh (zi) - [^oi + ^02 cosh {h) + sinh (^i)] cosh (zi)} ) , 
a V A2 



(23) 



where 



A2 = ((5oi -I- ^02) cosh [h) + (1 + 6oi5o2) sinh [h) , 



and mean value of susceptibility x'- 



1 



X = 



2 (cosh (/i) + {5m + ^02) sinh {W) 



{5oi + (^02) cosh (/i) + + 5m5o2) sinh {h 



,T>Tc. 



(24) 



The thickness dependence of mean susceptibility 
(Eq.(p4|)) is reported in Fig. 10. It is seen from compari- 
son of the Fig. 10 and left branches of the solid curves in 
Fig. 8 that x behaviour in paraelectric phase at T > Tc 
and at Td < T < Tc is strongly different: there is nei- 
ther divergency nor maximum at critical thickness, but 



X increases slowly with thickness increase. As a matter 
of fact this difference is related to the Ic temperature de- 
pendence: /c = at T > Tc- It is worth stressing that 
the considered case T > Tc corresponds to incipient fer- 
roelectric films, which are known to be in paraelectric 
phase with Tc = 0, so that Td = too {Td < Tc by 
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definition). 



2 



n(z) 

Ho 



0.6 




spends to I Ic or T ^ Td . 

The calculations can be performed by expansion of ex- 
pressions (14d) and ( |l9a| ) for fl and x respectively into 
power series in m up to the first nonvanishing terms. 
Because of complex form of n(m) arid y(m) (secAppen- 
dices 1, 2 and expressions (15a, 15b , 15c and |l9b| )) we 
put the detailed calculations in Appendix 3. The critical 
behaviour of dielectric susceptibility while approaching 
the transition fro m p araelectric phase is obtained directly 
from expression (22) that is valid for T^i <T < T^- These 
calculations lead to the following results: 

Mean pyroelectric coefficient: 



ni 



- i>i 



Vh — he 



(25a) 



z//o 

FI©i 4. Pyroelectric coefficient profiles at Soi — 0.5, 
6o2 = 2 for the following m values: 0.04; 0.12; 0.33; 0.57; 



0.99; 1 - 10" 
5: 6; 7. 



1 — 10 respectively for the curves 1; 2; 3; 4; 




FIG. 5. The thickness dependence of inverse mean pyro- 
electr'fc coefficient at T < Td for the following ((5oi, ^02) val- 
ues: (1.5, 2), (0.6, 0.6) and (0.03, 0.02) for solid curves 1, 2 
and 3 respectively. The dashed lines illustrate the description 
of inverse mean pyroelectric coefficient by the approximate 
expression (25a) with the same parameters as above. Arrows 
show the values of critical thickness. 



VI. CRITICAL BEHAVIOUR OF 
PYROELECTRIC COEFFICIENT AND 
DIELECTRIC SUSCEPTIBILITY 

In this section we consider the behaviour of pyroelec- 
tric coefficient and dielectric susceptibility in the vicinity 
of the thickness induced ferroelectric phase transition in a 
film. This means that we have to calculate the behaviour 
of these quantities at infinitesimal electric polarization, 
i.e. parameter m 0. On the other hand m corre- 



ni 



c 



py 



,T<r, 



cl ■) 



(25b) 




T/Te 

FIGp|g§. The temperature dependence of the inverse mean 
pyroelectric coefficient for the films with following dimention- 
less thickness A > A^: 2, 3 and 8 (see respectively solid lines 1, 
2 and 3) at di = 0.5, d2 = 2. The dashed lines illustrate the 
description of the inverse mean pyro electric coefficient given 
by the approximate expression ( 25b ) with the same parame- 
ters as above. The dash-dotted line illustrate the dependence 
of the inverse pyroelectric coefficient of thick film HoiT). Ar- 
rows show the values of critical temperature. 

Thickness dependence of mean dielectric susceptibility 
in paraelectric phase: 



-I 2Ci^ 

Ml^l^-O — 1 1 
tic — '1 



<Tc,l <lc, 



and in ferroelectric phase: 

<1 — tic 



,T<TciA>lc 



(26a) 



(26b) 



Temperature dependence of mean dielectric susceptibil- 
ity in ferroelectric phase: 



8 



x\ 



Tn-T 



cl ; 



and in paraelectric phase: 



(27a) 



(27b) 



where constants Cpy, Cpy, C^i^, can be found in Ap- 
pendix 3. 



15 
10 

5 

%(z) 

2 




2 4 6 8 

z/lo 

FIG. fZ? Dielectric susceptibility profiles in ferroelectric 
phase for the same parameters as those in Fig. 4. 




FIG. f8.s The thickness dependence of inverse mean dielec- 
tric susceptibility at T < Tc for the same parameters as those 
in Fig. 5 (solid lines 1, 2 and 3). The dashed lines illustrate 
the description of inverse mean dielectric susceptibility given 
by the approximate expressions ( 26a 261: ) . Arrows show the 
values of critical thickness. 



The dependencies given by expressions (25a- 27b| ) are 
shown by dashed lines in Figs. 5, 6, 8 and by crosses in 
Fig. 9 . It is seen that for thin enough films the expressions 
27a| , 27b) for the temperature critical behaviour, 
Tci) describe pretty good 



(& 



in the vicinity of t 
the temperature dependence of x and 11 far enough from 
these points, while for the thickness critical behaviour 



the approximate description ( p5a| , 26a, 26b) is good only 
in the close vicinity of Ic (see Figs. 5, 8). Note that the 
lines, corresponding to exact and approximate solutions 
for curves 1 in Figs. 9 and 6 are almost the same, i.e. ap- 
proximate formulas perfectly fitted 11 and x for the most 
thin films. This is because thickness induced ferroelectric 
phase transition is typical feature of thin films. 



1 



Xq(t=Q ) 

%(T) 
0.5 




0.5 



1 . 



1 . 5 



T/Te 



FIG. The temperature dependence of the inverse mean 
dielectric susceptibility for the films with for the same param- 
eters as those in Fig. 6 (solid lines 1, 2 and 3). The crosses 
illustrate the dependence of inverse mean dielectric suscepti- 
bility given by the approximate expressions (29a, b). Arrows 
show the values of critical temperature. 

It is important to emphasize tha t cri tical dependence 
of pyroelectric coe fficient ( sec Eq.( t25b| )) and dielectric 
susceptibility (Eg. ( p7a[|27b| ) ) on temperature is similar 
to that of bulk ferroelectrics. 

The critical behaviour of the films with different thick- 
ness is related to size effect of transition temperature 
Tci- This effect is described by e xpression (13a) and is 
reported on Fig. 3. Formula (13a) can be simp ified for 



thick enough film and for thin films with thickness about 
critical one, namely 



Tci 



A > 1 (28a) 




A - Ac < 1 (28b) 



It is seen from Eqs.( p8b| ) that at A — )■ oo, i.e. I — > oo, 
Tci 1, he. Tci ~> Tc, and at A — > Ac Td — > in agree- 
ment with behaviour shown in Fig. 3. It follows also from 
Fig. 3 that thick film limit can be obtained for finite A val- 
ues even for small enough di_2 (see curve 3, where Td w 1 
at A « 15 and di = 0, 2 , d2 = 0, 1). To show how the ap- 
proximate expressions ( ^8b| ) fit the exact curve (13a) we 
show them in Fig. 3 (for the parameters of curve 2 from 
that figure) by dotted and dashed lines respecti vely One 
can see the satisfactory agreement between Eqs.(^8b|) and 
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the curve 2, the fitting by expres sion (28b) is better (for 
all values of A) than that by Eg. (28a). 



VII. DISCUSSION. COMPARISON WITH 
EXPERIMENT 

We performed the analytical calculations of pyro- 
electric coefficient and dielectric response profiles, their 
thickness and temperature dependence on the base of 
solution of differential Lame equation. The considered 
boundary conditions include different values of extrapo- 
lation lengths on the film boundaries. To our mind this 
is more adequate for single film (not multilayer structire) 
than the supposition used in the previous works, where 
the same extrapolation length (i.e. the same conditions 
on both boundaries) have been used. In addition, the 
majority of previous works (see e.g. |^-[ic[|) have been de- 
voted to numerical calculations of the properties of thin 
films made from BaTiOa and PbTiOs materials. The 
analytical expressions for description of size effects of py- 
roelectric coefficient, dielectric response and polarization 
in ferroelectric thin films have the advantage of their ap- 
plicabilitv to anv film one^^can be interested in. 



. 8 



0.4 




10 



15 



Wo 



FIG. 10. The calculated thickness dependence of mean di- 
electric susceptibility at T > Tc for the films with the same 
parameters as those in Fig. 8. The inset shows the observed 
thickness dependence of real part of permittivity (measured 
at 100 kHz) of SrTiOs film @. 

The theory predicts the divergency of dielectric sus- 
ceptibility and pyrocoefficient in the vicinity of thickness 
induced ferroelectric phase transition. Temperature de- 
pendence of these quantities can be described by Curie- 
Weiss law in wide enough temperature range while in 
thickness domain the same dependencies x ~ 1/ K ~ 'c| 
and n ^ 1/ — Ic are valid only in the vicinity of crit- 
ical thickness Ic- To our mind these simple laws can be 
in real help for those interested in these properties mea- 
surements. 

Unfortunately the results of measurements of size ef- 
fects in ferroelectric thin films are strongly restricted. 



This may be related to the difficulties of thin films fabri- 
cation (especially with the thickness in the vicinity of 1^ 
that can be several tens nanometers), to the influence of 
electrodes on the film properties etc. (see e.g. |]l6|). 
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FIG. 11. Py r ocoe fficient profile calculated on the base of 

0.6, 



13c 



14a) in section 5 at m = 0.1; (5o 



formulas" 

<5o2 ~ 0.05 (solid line) and measured in PZT film |l^ (dashed 
line). 

To overcome the latter difficulties the application of 
spectroscopic methods and optical and radiospectroscopy 
methods in particular seems to be desirable. Special 
methods are needed also for investigations of the film 
properties profiles. The most informative is ellipsome- 
try method (see | p7| and ref. therein), making possible 
to measure optic refraction index profile. The calcula- 
tions of polarization profiles similar to those in Section 
3 made it possible to show that measured refraction in- 
dex profile is related to the film polarization profile . 
The special laser-intensity-modulation method (LIMM) 
can be applied to pyrocoefficient profile measurements 
(see [^,^). Pyrocurrent coefficient profile extracted in 
pSj from the pyrocurrent measurement by LIMM in 1 
/im thick PbZro.75Tio.25O3 film is shown by dashed line 
in Fig.ll. It is seen that it is looks really like some 
of the theoretical curves in Fig. 4. Detailed comparison 
of the theory with experiment is shown in Fig.ll where 
solid line is theoretical curve with parameters m = 0.1, 
di = 0.6, d2 = 0.05. We choose to fit the right "tail" of 
theoretical curve because of the absence of data related 
to high frequency measurements so that it can be ambi- 
guities in the experimental points near left " tail" bound- 
aries. However the procedure of pyrocoefficient profile 
extraction from the measured (in broad frequency range) 
pyrocurrent implies some numerical methods, the search 
for the most effective one is still under way (see e.g. [Q). 
For solution of this problem (that is the essential part of 
LIMM application to thin films) the calculations of 11 (z) 
profile seem to be important. On the other hand, spa- 
tial charge is present at the near-electrode layer in real 
ferroelectic films, that leads to the suppression of the fer- 
roelectic long range order in this region Ml. This phe- 
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nomenon has not been taken into account in the present 
work. Therefore both experimental and theoretical ap- 
proaches need further improvement for the better com- 
parison of the theory with the experiment. 

Unfortunately there is no measurements of pyroelec- 
tric coefficient thickness dependence. Static dielectric 
susceptibility thickness dependence measured recently in 
SrTiOa thin films [|6|; it is shown on the inset to Fig. 10. 
Since SrTiOs is an incipient ferroelectric, i.e. it is in 
paraelectric phase for all temperatures {T > 0), x be- 
haviour is described by Eq. (p3) (compare curves in the 
Fig. 10 with those in the inset jpretty good. 

The theory prediction about the divergency of pyro- 
electric coefficient and static dielectric susceptibility in 
the vicinity of thickness induced ferroelectric phase tran- 
sition in thin films is waiting for the experimental confir- 
mation. 



1 + - m{l + m) yo(x) 

1 — TO 



(A4) 



To calculate the integral in the expression ( A3d ) it is 
necessary to know the explicit form of the function gy (x) . 
This function has the form (see Eqs.(^, pa|)): 



gpi{x) = -2(1 + TO)nr 



2to 



(A5a) 



1 + 

for the calculation of pyroelectric coefficient 11(2;) and 

5,,,(a;) = -2(1-1- to)xo (A5b) 

for the calculation of dielectric susceptibility x{x)- After 
substitution (A_5a,A5b) to (A3d) and integration over x 
and X 1 21 1, one can get the following expressions for the 
inhomogeneous solution of equation ( [A1| ) : 



APPENDIX A: 



y^ix) - 



^ 2 i?(am(.T), to) \ , , 
X yi{x) 

, 1 — TO ' 



Let us rewrite the equations (^^ and ( |5a| ) for dielec- 
tric susceptibility x{x) and pvroelectric coefficient Il{x) 
with respect to expression ^ for polarization distribu- 
tion. Introducing new variable a; = {zi + zq)/^/! + m and 
denoting x ^-^d 11 as one can find following equation 
for y{x): 



<Py{x) 
dx^ 



+ {l + m-6myQ{x))y{x) = gy{x), (Al) 



with the boundary conditions 
dy{x) 



dx 
dy{x) 



yo{x) = sn(a;,TO) 



dx 



302 



y{x)\x=x 



(A2) 



General solution of the equation ( Alj) (particular case of 
Lame equation) , is the following [ 20[ : 



y{x) = Ci yi{x) + C2 y2{x) H- y^ix) 
yi{x) — cii{x, m) dii{x, m) 
dx 



y2{x) = yi{x) 
y2{x) = yi{x) 



{yi{x)f 



9v {x) yi jx) 

{yi{i)f 



dxdx 



(A3a) 
(A3b) 

(A3c) 
(A3d) 



The integral in the expression ( A3c) can be easily reduced 
into a standa rd fo rm (see |2^ ) , its further substitution to 
the formula (A3c) gives the explicit form of the second 
fundamental solution of the equation ( |Al[ ): 



y2{x) 



-E{am{x),m) ) yi{x 



+2/0 (a;) 



1 + TO — 2to 1/9(2;) 



1 — m 



n, 



^2to(1 +to) 

1 — TO 



(A6a) 



and 



rh , ^ X 1 + m — 2m y?:(x) , , „, ^ 

ytix)^~2{l + m)xo (A6b) 

(1 — to) 

Constants Ci , C2 are determined from the bou ndary 
conditions (A2). Substituting expression ( A3a ) into 
( [A2| ), one can obtain the system of linear equations for 
calculation of Ci and C2. The solution of this system 
gives following form for constants Ci and C2 



Ci = 



c. = 



BIA22 — B2A21 
A ' 

B2A11 — B1A12 
A ' 



A = ^22^11 - A21A1 



(A7a) 



where 



All = 
A21 = 
A12 = 
A22 = 



VI + 


TO 




VI + 


TO 


Soi 


VI + 


TO 


Sqi 


VI + 


TO 



Soi 



yiixi) - y[{xi), 
y2{xi) - y2{xi), 
yi{x2) + y[{x2), 
2/2(2^2) + 2/2(3^2), 



Bi = y3[xi) J 2/3(2^1), 



B2 = -2/3(2^2) 



5oi 



vr 



301 



-2/3 (a;2) 



(A7b) 



11 



Using the explicit form of solutions yi{x), y2{x), ysix) 
( A3b , A4 , A6a , A6t), values of xi, X2 and notations for 
ei and 62 

xi — F(arcsin(/i), m), X2 ~ 2K{m) — F(arcsin(/2), m) 
ei — i?(arcsin(/i), to), 62 ~ 2E{m) — i?(arcsin(/2) ^ 



,TO 



1 + TO 

2to 



4m 



(1 + to)^ 



it is easy to obtain following expression for Aij and Bi 

All = /i ("(1 + "i) (1 + ^) - 2™ A") , 



^12 = -/2 ((1 + to) (^1 + i-^ - 2m /I 
1 



A 



2i 



1 — m 



^ VI + TO 2m + m(l + m)ff 
' Soi 1 - TO 



, 1 + TO 

-Au X, - e 

1 — TO 



(A8) 



So 
2e,; 

1 — TO 



1 — TO 



B: 



tch 



i = 1,2. 



2(1 + to) Vl + TO. , „ 

^ ^ (1 + m + 2to // 

(1 - to) (5oi 



(A9) 



(AlO) 



This permits to get final for m of coeffi cien ts Cf, Cf*' 
with the help of expressions (A7a, |A8| , A9 ) a nd c oeffi- 
cients Cf*', with the help of expressions (A7a, AS, 



aTo| ) 



APPENDIX B: 

Mean pyroelectric coefficient 11 and dielectric suscepti- 
bility X E^i'e introduced in the expression (^. After chang- 
ing of variables in the integrals (^), the mean value of 
function g(x) can be written as 



the exphcit form ( A3b| , A4 , A6a , A6t) of yi(x), y2{x), 
2/3 (x), it is easy to reduce the integrals to the standard 
form and to find that (see |2^]) 



y2{x)dx 
1 

1 — TO 

yl^{x)dx 



yi{x)dx = 2/0 (a;). 



1 + m \ 1 + m 
X- E(am{x),m) yQ(x) - yi{x) 



1 — m 



1 — to' 



2 i?(am(a;), m) \ 2yi(x) 
X z )yo{x) 



1 — m 



1 — m 

ch/ X , r, 1 + 2£'(am(a;),m)\ 

yT{x)dx = 2x0:^ [x -r-^ — - ■ 

1 — TO\ 1— m / 



These expressions along with values of a;i and X2 per- 
mit to get the mean values of yi{x), 2/2(2;), 2/3 (x) in the 
following form 



Vi 



/2-/1 

X2- Xi 



(B2a) 



1 



y2 = 



{x2 - xi) (1 - m) 



. , 1 + m / VT + m 

J2 a;2 - -; I 62 



1 — m 



302 



—PV 

Vs = 



1 — TO 

n. 



1 + TO / vT+m 

-JlXl - I 61 ' 



301 



(B2b) 



{X2 - Xi) (1 - to) 



/2(2(e2-^^^^)-(l-m)x2)- 



-A(2(6i-^^li±^)-(l-m)xi 



yc^^^^l±^(^,.^, 2(62-61) 



a;2 — xi 1 — m 



1 — m 



(B2c) 
(B2d) 



X2 - Xi 



g{x) dx 



Using expressions ( |13c| , |17[ ) for x(x) and n(x), we can 
write their mean values as follows 



n^Ciyy, + Ciyy2 + yf 

x = cf'm + cty2 + yt 



(Bla) 
(Bib) 



To calculate 11 and x it is necessary to integrate all the 
solutions of equation (Al) over x. Taking into account 



1. Ferroelectric phase 

For the investigation of behavior of mean pyroelectric 
coefficient fl and dielectric susceptibility x in the vicinity 
of thickness induced phase transition, we are to expand 
n(TO) and y(m) at to — > 0. Using the series for elliptic in- 
tegrals [ pl[, on e can expand Ay , Bi and (see (A8, A£, 
A10| ) and ( |B2a| , [A34 |B24 |B2d| )) to the first nonvanishing 
terms and obtain following relations 
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n 



\Soi 



S2 

^05 



where 



3Ao/ 



Ic 



(B3a) 



(B3b) 



lie = TT — arcsin (Si) — arcsin (6*2) 



Ao 



The dependence of dimensionless thickness li{m) has the 
form 



h — \/l + m {x2 — xi) . 
The expansion over m gives following relation 



kL^a = he + ^Aom 



(B4) 



Substitution of this expression into ( B3a , B3b ) gives the 
structure of mean pyroelectric coefficient 11 and dielectric 
susceptibility x i^i the vicinity of thickness induced phase 
transition — > + : 



^ch 



XUi^hc+O 7 _ 7 

il Ll 



\Soi 



S2 
<5o; 



he 



n 



CL 



CL = n 



+ 



V3A^/i 



(B5a) 
(B5b) 

(B6a) 
(B6b) 



The temperature dependence of 11 and x in the vicin- 
ity of the phase transition can also be derived from (B4) 
and (B3a, B3b) with respect to obvious form of li and 
Soi- Introducing new notations 



h = AVT^,A = ;i(r = 0), 

Soi = diVl - r, di = (5m (T = 0), 
A, = ;ie(T = 0),T = |^, 



(B7) 



and using equation for critical temperature Td = T^/Tc 



A\/l — ''c/ = TT — arcsin 



v/l + d?(l-r,z). 



— arcsm 



Vl + di(l-Tez) 



we can expand the expression ( B'l ) and get following se- 
ries up to the first order in m 



1(1- T,i) Ao (r = T,i) 



\ /-I _ I di yi — Tei I d2\/l-Tci 



Applying this expansion to (B3a, B31:), we can rewrite it 
as follows 



XIt-»t„,--0 



Tel - T 

2N^ 



XoiT = 0), 



n 



c 



PV 



(B8a) 
(B8b) 

(B9a) 
(B9b) 



where 



N 



= AcVl - 



di^/1 - Tcl 



0^2 \/l - Tc, 



l + rff(l-T,0 l + rf2(l_r,0' 



3 (1 + rf2 (1 _ ^^,))^ 3(l+d2(i_^^^))^ 



2. Paraelectric phase 

Mean dielectric susceptibility x in paraelectric phase 
(i.e. at temperature above Tc or at h < he) has the form 

_ ^ /I (^oi+^o2)sin(/i) + 2(l-cos(^i)) _ ^ 
° Vi (1 - 6016Q2) sin (?i) + (Soi + S02) cos (/i) 

(BIO) 



Expanding (BIO) in terms of (h ~ he) one can write the 
asymptotic form of x at /i ^ he 



2Ci 



ch 



he — h 



(Bll) 
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where C^^ is given by expression ( B5b| ) . Using the tem- 
perature dependence (B7) of li and 6oi, it is easy to 
rewrite (BIO) as follows 



T -Tel' 



(B12) 



where coefficient is given by the expression ( B8b|) 
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